Interferometric determination of the s- and d-wave scattering amplitudes in ^''Rb 
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We demonstrate an interference method to determine the low-energy elastic scattering amplitudes 
of a quantum gas. We linearly accelerate two ultracold atomic clouds up to energies of 1.2 mK and 
observe the collision halo by direct imaging in free space. From the interference between s- and d- 
partial waves in the differential scattering pattern we extract the corresponding phase shifts. The 
method does not require knowledge of the atomic density. This allows us to infer accurate values 
for the s- and d-wave scattering amplitudes from the zero-energy limit up to the first Ramsauer 
minimum using only the Van der Waals Ce coefficient as theoretical input. For the *^Rb triplet 
potential, the method reproduces the scattering length with an accuracy of 6%. 

PACS numbers: 34.50.-s, 32.80.Pj, 03.75.-b, 03.65.Sq 
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The scattering length a, the elastic scattering ampli- 
tude in the zero-energy limit, is the central parameter in 
the theoretical description of quantum gases Q, Q, y| . It 
determines the kinetic properties of these gases as well as 
the bosonic mean field. Its sign is decisive for the collec- 
tive stability of the Bose-Einstein condensed state. Near 
scattering resonances, pairing behavior 'z| and three- 
body lifetime |3 can also be expressed in terms of a. 
As a consequence, the determination of the low-energy 
elastic scattering properties is a key issue to be settled 
prior to further investigation of any new quantum gas. 

Over the past decade the crucial importance of the 
scattering length has stimulated important advances in 
coUisional physics l^]. In all cases except hydrogen |3] 
the scattering length has to be determined experimen- 
tally as accurate ah initio calculations are not possible. 
An estimate of the modulus \a\ can be obtained relatively 
simply by measuring kinetic relaxation times l6l| . In some 
cases the sign of a can be determined by such a method, 
provided p-wave or d-wave scattering can be neglected 
or accounted for theoretically ^7]. These methods have a 
limited accuracy since they rely on the knowledge of the 
atomic density and kinetic properties. Precision deter- 
minations are based on photo-association ^] , vibrational- 
Raman |3| and Feshbach-resonance spectroscopy |ldllll |. 
or a combination of those. They require refined knowl- 
edge of the molecular structure in ground and excited 
electronic states |j| . 

In this Letter we present a stand-alone interference 
method for the accurate determination of the full (i.e. 
complex) s- and d-wave scattering amplitudes in a quan- 
tum gas. Colliding two ultracold atomic clouds we ob- 
serve the scattering halo in the rest frame of the coUi- 
sional center of mass by absorption imaging. The clouds 
are accelerated up to energies at which the scattering pat- 
tern shows the interference between the s- and d- partial 
waves. After a computerized tomography transforma- 
tion J12l | of the images we obtain an angular distribu- 
tion directly proportional to the differential cross sec- 



tion. This allows us to measure the asymptotic phase 
shifts rii{k) (with k the relative momentum) of the s- 
wave (/ — 0) and d-wave [l = 2) scattering channels. Us- 
ing these ?7;(fc) as boundary conditions, we integrate the 
radial Schrodinger equation inwards over the —Cq/t^ tail 
of the potential and compute the accumulated phase [l3| 
of the wavefunction at radius 20 oq (with ap the Bohr 
radius). All data of •qi{k) are used to obtain a single op- 
timized accumulated phase from which we can infer all 
the low-energy scattering properties, by integrating again 
the same Schrodinger equation outwards. Note that this 
procedure does not require knowledge of the density of 
the colliding and scattered clouds, unlike the stimulated 
raman detection approach of Ref. [lj| . We demonstrate 
this method with *^Rb atoms interacting through the 
ground-state triplet potential. We took data with both 
condensates and thermal clouds. Here we report on the 
condensates, as they allow to observe the largest range 
of scattering angles, 25° < 6* < 90°. Up to 80% of the 
atoms are scattered without destroying the interference 
pattern. With our method, we obtain a = -1-102(6) ao for 
the scattering length. The d-wave resonance [lj| is found 
at the energy Eres = 300(70) /xK. These results coincide 
within experimental error with the precision determina- 
tions (a = 98.99(2) ao [HEl and Eres = 270 /iK [l^), 
obtained by combining the results of several experiments 
as input for state-of-the-art theory. 

In our experiments, we load about one billion *^Rb 
atoms in the (fully stretched) |i^ = 2, tof = 2) hyper- 
fine level of the electronic ground state from a magneto- 
optical trap (MOT) into a loffe-Pritchard quadrupole 
trap (21 X 477 Hz) with an offset field of Bq = +0.9 G. 
We pre-cool the sample to about 6 /xK using forced radio- 
frequency (RF) evaporation. The cloud is split in two 
by applying a rotating magnetic field and ramping Bq 
down to a negative value Bq . This results in two Time- 
averaged Orbiting Potential (TOP) traps loaded with 
atoms [ijj. By RF-evaporative cooling we reach Bose- 
Einstein condensation with about 10^ atoms in each 
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density distribution of the halo in cylindrical coordinates, 
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FIG. 1: a) Optical density of the scattering halo of two Rb 
condensates for collision energy E/ks = 138(4) /jK, measured 
2.4 ms after the collision (upper half: measured; lower half: 
calculated lig|| ); b) radial density distribution obtained af- 
ter tomography transformation of image-a (upper half: mea- 
sured; lower half: calculated [lyDi c) the dots show the angu- 
lar scattering distribution W (0) obtained after binning plot 
b, the line is the best parabolic fit. d,e,f) As in plots a,b,c, 
but measured 0.5 ms after a collision at 1230(40) /jK. The 
field of view of the images is ~ 1 x 1mm . 



cloud and a condensate fraction of ~ 60%. 

We then switch off the TOP fields and ramp Bq back 
to positive values, thus accelerating the clouds until they 
collide with opposite horizontal momenta at the location 
of the trap center. The collision energies E = \2^bBq \ = 
h'^k'^/m (with /XB the Bohr magneton and m the mass of 
^"^Rb) range from 138 /xK to 1.23 mK with an overall un- 
certainty of 3% (RMS). Approximately 0.5 ms before the 
collision we switch off the trap. A few ms later we observe 
the scattering halo by absorption imaging. Fig.^ (upper 
part) displays the s- wave-dominated scattering halo (av- 
eraged over 20 pictures) of fully entangled pairs (see |l8l |) 
obtained for a collision energy of E/ks — 138(4) /iK. In 
Fig.ni (upper part), taken at E/ks = 1.23(4) mK the 
halo is entirely different, showing a d- wave-dominated 
pattern. The lower halves of Fig.^ and Fig.QJi show the 
theoretical column densities n2 (x, z) — J n{x, y, z) dy, 
where n (x, y, z) is the calculated [lS| density of the halo. 

As the atoms are scattered by a central field, the scat- 
tering pattern must be axially symmetric around the 
(horizontal) scatterin g ax is (z-axis). As pointed out by 
the Weizmann group |20j , this allows a computerized to- 
mography transformation tij to reconstruct the radial 



1 /? 
Here p — {x^ -I- y^) , n2 {nx^ z) is the ID Fourier trans- 
form along the a;-direction of the optical density with 
respect to z, and Jq {q) is the zero-order Bessel function. 
The transformed plots corresponding to the images of 
Fig.^,d are shown as Fig.n3,e respectively. 

To obtain the angular scattering distribution W{6) the 
tomography pictures are binned in 40 discrete angular 
sectors. For gas clouds much smaller than the diam- 
eter of the halo, W (6) is directly proportional to the 
differential cross section a (9) = 2tt \f (9) + f {tt — 6)\ . 
Here, the Bose-symmetrized scattering amplitude is given 
by a summation over the even partial waves, / (9) + 
fin -9)^ {2/k) Ei=e.en(2«+l)e""f/(cos^) sin 77,. Note 
that unlike in the total elastic cross section [a = 

J^^'a{0)sm9d9 = (STr/fc^) ^,^_J2; + l)sinSO, 
the interference between the partial waves is promi- 
nent in the differential cross section. Given the 
small collision energy in our experiments, only the s- 
and d-wave scattering amplitudes contribute, fs (9) + 
fs{n-e) = (2/fc)e^''«sinr;o and fd[9) + fdi^ - 9) = 
(2/A:)(5/2)e*''^ (3cos2 9-1) sin 7^3. Therefore the differ- 
ential cross section is given by 



o-(6') = — sm ?7o 



25 
1 + 5 cos(77o — rj2)u + —u^ 



, (2) 



where u = (sin 772 /sin 770) (3cos^ 6 — 1). 

To obtain the phase shifts, we plot the measured an- 
gular distribution W [9) as a function of (3 cos^ 9 — \) as 
suggested by Eq. 10). The results for Fig.^ and Fig.QJi 
are shown as the solid dots in Fig.^ and Fig.^, re- 
spectively. A parabolic fit to W [9) directly yields a 
pair {"q^^ik), 77™^(A:)) of asymptotic phase shifts (defined 
modulo tt) corresponding to the two partial waves in- 
volved 21]. The absolute value of W [9) depends on 
quantities that are hard to measure accurately (like the 
atom number) so we leave it out of consideration. We 
rather emphasize that the measurement of the phase 
shifts is a complete determination of the (complex) s- 
and d-wave scattering amplitudes at a given energy. 

The radial wavefunctions corresponding to scattering 
at different (low) collision energies and different (low) 
angular momenta should all be in phase at small in- 
teratomic distances [13. This so-called accumulated 
phase common to all low-energy wavefunctions can be 
extracted from the full data set {{r]^^{k),'q^^{k))} men- 
tioned above. In practice, we use the experimental 
phase shifts 77Q'^^(fc) and ri^^{k) as boundary condi- 
tions to integrate inwards - for given E and I - the 
Schrodinger equation h^ d^ xi''') I dr"^ +P^(^)x('') = 0, and 
obtain the radial wavefunctions xi'"')/''' down to radius 



nil 



w 

CD 



-71/2 



Till 



CD 
D. 



-ji/2 



a) d-wave 


i 


\ 




** 


\ 



4x10 



10 



100 



1000 



= 1 — 1 1 — 1 

b) s-wave 


\i 


\ 


"^^^^^ 


H 


' 


■ 



10 100 1000 Ert 

Collision Energy (p,K) 



FIG. 2: a) d-wave and b) s-wave phase shifts versus coUision 
energy in /iK. The circles are the results of the parabohc fit of 
W{6) for individual images. The full black lines is calculated 
from the accumulated phase "I>opt optimized from all data 
points. The grey lines show the influence of the uncertainty 
of ±7r X 0.025 on "l>opt. (The vertical dotted line indicates the 
condition ?7o = »?2). s — d interference is only observed in the 
gray areas. The first s-wave Ramsauer-Townsend minimum 
is found at Ert = 2.1(2) mK. 



ri„ = 20ao. Here, p2(r) = m{E -Vir)) - hH{l + l)/r'^ , 
where V{r) ~ —Cq/t^ is the tail of the interaction poten- 
tial. At radius 20 ao, the motion of the atoms is quasi- 
classical and the accumulated phase can be written as 
<l>(r) ~ arctan [p(r) /(^ d In ~x /dr)] . The distance 20 ao is 
small enough J23] for 4>(n„) to be highly insensitive to 
small variations in i? or I [l3j and large enough that the 
— Ce/r^ part of the interaction potential is dominant over 
the full range of integration. With a least-square method 
we establish the best value $opt('''m) = 1-34 ± tt x 0.025 
for the accumulated phase at 20 ao |23| . Here the error 
bar reflects the experimental accuracy and not the sys- 
tematic error related to the choice of Cg , the latter being 
of less relevance as discussed below. Interestingly, the 
d-wave scattering resonance [15J results in a sudden vari- 
ation of 7^2'^'' with the collision energy in the vicinity of 
that resonance (see Fig. |5^). This imposes a stringent 
condition on the optimization of ^opt and constrains its 
uncertainty. 

Once $opt has been established, one can use it as a 
boundary condition to integrate the Schrodinger equation 
outwards and compute 'qi{k) for any desired (low) value 
of k and /. Fig.|21 shows the resulting phase shifts for 
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FIG. 3: s-wave (dashed line), d-wave (dotted line) and total 
(full black line) elastic cross sections (in cm^) versus collision 
energy (in /^K), computed from the optimized accumulated 
phase $opt as determined in this work. The gray lines are the 
total elastic cross sections, obtained from $opt ± tt x 0.025. 



collision energies up to 5 mK |2J|. The first Ramsauer- 
Townsend minimum |25| in the s-wave cross section is 
found at coUision energy Enr/kB = 2.1(2) mK. The solid 
dots represent the rfi^^{ki) obtained from the parabolic fit 
of W [6) from individual images. The three open circles 
correspond to measurements for which the sign of the 
phase shifts could not be established ;26|. Refinements 
to the present data analysis may include the occurrence of 
multiple scattering as well as the influence of the spatial 
extension of the colliding clouds taking into account the 
non condensed fraction. 

Knowing the phase shifts, we can infer all the low- 
energy scattering properties. Our results for the elas- 
tic scattering cross section are shown in Fig.|3| The 
(asymmetric) d-wave resonance emerges pronouncedly 
at 300(70) /iK with an approximate width of 150 jjK 
(FWHM). Most importantly, the scattering length fol- 
lows from the A: — > limiting behavior, rj^ik — > 0) = —ka. 
We find a = +102(6) an, whereas the state-of-the-art 
value is a = 98.99(2) ao [ia |. 

Comparison with the precision determinations |lll llq 
shows that our method readily yields fairly accurate re- 
sults, relying only on input of the Cg coefficient. We 
used the value Ca = 4.698(4) x 10^ a.u. ^. In the 
present case, one does not need to know Cg to this accu- 
racy. Increasing Cg by 10% results in a 1 %-change of the 
scattering length. Clearly, the systematic error in <i>opt 
accumulated by integrating the Schrodinger equation in- 
ward with a wrong Cq largely cancels when integrating 
back outward. However, in the case of a s-wave resonance 
other atomic species may reveal a stronger influence of 
Cg on the calculated scattering length. Simple numerical 
simulations show that the value of Cg becomes critical 
only when the (virtual) least-bound state in the interac- 
tion potential has an extremely small (virtual) binding 
energy (less than 10~^ level spacing). Hence our method 



should remain accurate in almost any case. 

This method can therefore be applied to other bosonic 
or fermionic atomic species, provided the gases can be 
cooled and accelerated in such a way that the lowest- 
order partial-wave interference can be observed with 
good energy resolution. We speculate that the accuracy 
of the method can be strongly improved by turning to 
smaller optical-density clouds and fluorescence detection. 
It will enable higher collision energies and observation of 
higher-order partial-wave interference. The use of more 
dilute clouds and longer expansion times will also elim- 
inate multiple-scattering effects and finite-size convolu- 
tion broadening of the interference pattern. Finally it will 
enable precision measurements of the scattered fraction, 
which in the case of ^'^Rb will allow us to pinpoint the 
location of the d-wave resonance to an accuracy of 10 fiK 
or better. In combination with state-of-the-art theory 
such improvements are likely to turn our approach into 
a true precision method. 

Similar experiments were reported during the final 
stage of completion of this Letter [23| . 
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number of data points, and the ??°''*(fci) are the phase 
shifts computed from "^opt- 

The p-wave {I = 4) phase shift can also be computed 
from the same ^opt, but the error made by assuming a 
constant accumulated phase increases like l{l -\- 1), and 
the resulting g-wave would be accordingly less accurate. 
Ramsauer-Townsend minima are observed around colli- 
sion energies where a phase shift crosses zero (see Fig.|5J . 
Since Eq. Q is unchanged under reversal of the sign of 
both 770"^ (fc) and 7)2^'^ (k), the result of the parabolic flt 
is not unique. We eliminate this ambiguity by compar- 
ing the diff'erence in accumulated phases obtained for 
Tj^^ik) and r{^^{k) with that obtained for —77™^ (fc) and 
— 77™''(fc). We find that in almost all cases, this difference 
- which should be vanishingly small - is much larger for 
one choice of signs than for the other. We conclude that 
the correct sign for a pair (»?o^''(fe), f?2^'^(fc)) is the one for 
which the difference between the accumulated phases is 
the smallest. In all but three cases (at the same collision 
energy), this criterium is conclusive. These three mea- 
surements (the open circles in Fig.|2l are left out of the 
procedure used to compute $opt. In hindsight they turn 
out to correspond to the marginal case mentioned in |21|. 
N. R. T homas, N . Kjargaard, P. S. Julienne, and A. C. 
Wilson, |cond-mat/0405544| 



